Consistent histories for tunneling molecules subject to collisional decoherence 
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The decoherence of a two-state tunneling molecule, such as a chiral molecule or ammonia, due to 
collisions with a buffer gas is analyzed in terms of a succession of quantum states of the molecule 
satisfying the conditions for a consistent family of histories. With hw the separation in energy of 
the levels in the isolated molecule and 7 a decoherence rate proportional to the rate of collisions, 
we find for 7 > u (strong decoherence) a consistent family in which the molecule flips randomly 
back and forth between the left- and right-handed chiral states in a stationary Markov process. For 
7 < uj there is a family in which the molecule oscillates continuously between the different chiral 
states, but with occasional random changes of phase, at a frequency that goes to zero at a phase 
transition 7 = to. This transition is similar to the behavior of the inversion frequency of ammonia 
with increasing pressure, but will be difficult to observe in chiral molecules such as D2S2. There are 
additional consistent families both for 7 > ui and for 7 < uj. In addition we relate the speed with 
which chiral information is transferred to the environment to the rate of decrease of complementary 
types of information (e.g., parity information) remaining in the molecule itself. 



I. INTRODUCTION 

The decoherence produced by the interaction of 
a quantum system with its environment is ubiqui- 
tous in nature and plays an important role in cur- 
rent quantum theory in at least two ways. First, 
it is widely believed that decoherence helps under- 
stand how the classical physics of macroscopic ob- 
jects emerges as an approximation to underlying 
quantum mechanical laws. Second, decoherence is 
the great enemy of quantum computation, quantum 
cryptography, and other schemes seeking to utilize 
specifically quantum effects for particular processes. 
For both reasons it is important to study specific mi- 
croscopic models from which one can hope to obtain 
general principles for decoherence. The present pa- 
per is the study of a simple two-level system which 
can be thought of as a crude microscopic model of 
chiral molecules or ammonia in which the lowest 
quantum energy levels correspond to the nearly de- 
generate eigenstates of a double- well potential, with 
decoherence occuring through collisions with parti- 
cles in the environment. 

Microscopic studies of decoherence are often 
framed in terms of a master equation for the den- 
sity operator of the decohering system. Such de- 
scriptions are perfectly valid, but because they rep- 
resent the average of a large ensemble of nominally 
identical systems, each with a different specific time 
development, they provide less information and less 
physical insight than the actual history of a single 
system. For example, in the phenomenon of inter- 
mittent fluorescence a single ion in a trap shows in- 
termittent light and dark periods when it does or 
docs not scatter resonance radiation This be- 
havior is not directly reflected in the density oper- 



ator, even though from the latter one can deduce 
parameters which govern the statistical behavior of 
the individual ion. 

Another way to understand the limitations of the 
density operator description is to consider its clas- 
sical analog for a Brownian particle confined to a 
small but macroscopic volume of a fluid by rigid 
walls. The probability distribution p(r, t) of the par- 
ticle position r will eventually tend to a constant 
over the region accessible to the particle, whereas 
the particle itself will continue to exhibit a sort of 
random walk. More details of what is going on in 
this steady-state situation is provided by the joint 
probability distribution of the sequence of successive 
positions n, 1*2, . . . of the particle at a sequence of 
times ti, ti, . . ., that is, its history. Averaging over 
a large number of histories will yield p(r,t), but in 
the process the information needed for a more de- 
tailed temporal description of the particle is lost. In 
the quantum case unravelings of the master equation 
provide a more detailed description of the micro- 
scopic time development, but these are often viewed 
as mathematical artifacts having no necessary con- 
nection with what is really going on in the quan- 
tum system. There are many possible unravelings; 
which, if any are correct? Standing in the way of an- 
swering this question is the infamous measurement 
problem of quantum foundations: textbook quan- 
tum mechanics introduces probabilities by means of 
measurements, but cannot say what it actually is 
that is being measured. 

However, the consistent histories or decoherent 
histories — hereafter simply referred to as histories — 
formulation of quantum mechanics, has no measure- 
ment problem, and provides the tools needed to 
identify trajectories or sequences of events that actu- 
ally correspond to physical processes. Or, putting it 
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another way, it allows one to identify certain classes 
of microscopic stochastic processes which can be 
consistently described in a fully quantum mechanical 
terms. The histories approach has previously been 
applied to quantum optical systems by Brun 0, H[ > 
though we believe the material presented here is the 
first application to the case of tunneling molecules, 
including chiral molecules. 

Early in the development of quantum mechanics 
the question arose as to why chiral molecules are ob- 
served in left- and right-handed versions even though 
the quantum ground state should be a symmetrical 
combination of the two forms. Hund Q provided the 
first step in addressing this paradox when he pointed 
out that the two enantiomers correspond to the two 
wells of a symmetrical potential with two minima, 
and that the time required to tunnel from one well 
to another for a typical chiral molecule is extremely 
long. A second step was provided by Simonius |5| 
who observed that interaction with the environment 
of a suitable sort (i.e., decoherence, though when 
he wrote that term was not yet current) can stabi- 
lize the chiral states for periods substantially longer 
than the tunneling time. At present it seems widely 
accepted that such decoherence is an important as- 
pect of the stability of chiral molecules, though there 
have been dissenting voices, e.g. @. 

The time dependence of the two-state model in- 
troduced in Sec. HH when analyzed in terms of con- 
sistent histories using the principles discussed in 
Sec. IIII1 and applied to specific consistent families 
in Sec. IIV1 yields some insight into this stability 
problem. In particular, we find that if the rate of 
decoherence due to collisions 7 (a parameter in our 
model) is much larger than the tunneling rate u in 
an isolated molecule, there is a consistent family in 
which the molecule spends a long but random pe- 
riod of time in each of the chiral states before flip- 
ping to the one of opposite chirality, in a two-state 
Markov process. As 7 decreases the flips become 
more rapid and the "dressed" quantum states be- 
tween which the flips occur become less and less chi- 
ral, with this type of family finally disappearing at 
a phase transition 7 = w. For 7 < to there is a dif- 
ferent consistent family with a rapid but continuous 
oscillation of the molecule back and forth between 
its chiral states, interrupted at random times by a 
change in phase. There are a variety of other con- 
sistent families, and these are discussed, along with 
their physical interpretation, in Sec. IIVI Most chiral 
molecules in most circumstances will be in the strong 
decoherence regime. We give some approximate nu- 
merical values in Sec. I VI Al for D2S2 in a buffer gas 
of helium, as it has been the subject of some careful 
calculations in Q. On the other hand the ammonia 
molecule, which though not itself a chiral molecule 
can behave like one in certain rotational states, has 
an inversion (tunneling) transition with a frequency 



that goes to zero with increasing pressure. This is 
probably an example of, or at least very similar to, 
the 7 = to phase transition, for reasons discussed in 
Sec. EH] 

An intuitive way of characterizing the stochas- 
tic time development in the histories formalism is 
to quantify the information dynamics, e.g., the loss 
over time of information about the system's ini- 
tial state due to random hopping. Such informa- 
tion loss is ultimately connected to decoherence [||, 
which can alternatively be viewed as flow of infor- 
mation about the system to the environment 0, 
and in Sec. [V] we illustrate the quantitative connec- 
tion between these two views of decoherence for our 
model. In our model decoherence corresponds to a 
flow of chiral information — i.e., is the molecule left 
or right handed? — to the environment. In Sec.lVlwc 
analyze this using quantitative measures defined in 
Scc. lIII Cl and compare the flow of chiral information 
to the environment with the decrease of complemen- 
tary types of information (e.g., parity information) 
about the earlier state of the molecule that remain 
in the molecule itself at later times. 

Our conclusions are summarized in Sec. I VIII 
which also indicates some ways in which the results 
reported here could be usefully extended. A few 
mathematical derivations and details are placed in 
appendices. 



II. MICROSCOPIC MODEL AND MASTER 
EQUATION 

A. Double-well potential and collisions 

We consider a quantum system, the molecule, 
with a double-well potential in which the two low- 
est energy eigenstates, |0) (even parity) and |1) (odd 
parity) , are sufficiently well separated in energy from 
all the higher levels that the latter can be ignored. 
The Hamiltonian is of the form 

H = (l/2)hwZ, (I) 

where Z = |0)(0| — is the Pauli operator 03 = 

a Zl so the energy splitting between the levels is fku. 
The linear combinations 

|i?) = (|0) + |l))/V2, |i) = (|0)-|l»/V2, (2) 

represent the right- and left-handed chiral forms of 
the molecule, or in ammonia the nitrogen on one or 
the other side of the plane formed by the hydrogens. 
In real molecules there are, of course, additional de- 
grees of freedom — rotations, vibrations, etc. We are 
assuming that for our purposes these can be ignored, 
i.e., the Hilbcrt space can be approximated as a ten- 
sor product of these other degrees of freedom with 
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the two levels representing the tunneling, with neg- 
ligible coupling between them. Hence the isolated 
molecule can be thought of as oscillating or tunnel- 
ing between the \R) and \L) states at an angular 
frequency w. In the Bloch sphere picture the kets 
\R) and \L) correspond to the points on the positive 
and negative x axis, and the sphere rotates about 
the z axis as time increases. 

Next we assume the molecule collides randomly 
with other particles (atoms or molecules), and the 
duration of each collision is short compared to the 
times we are interested in. Successive collisions need 
not be independent of each other, but we assume 
that correlations die away rapidly after some cor- 
relation time t c , which could be shorter than the 
average time between collisions in a dilute gas, but 
might be significantly longer in a dense gas or liq- 
uid. We will consider properties of the molecule at a 
succession of times to, t\, ti. . . , where the m'th time 
interval, A£ m = t m +i — t m is always greater than r c , 
and ideally should be significantly greater than r c . 
That is, we are using a description which is coarse- 
grained in time; the importance of this will appear 
later. During the m'th time interval there may be 
zero or one or more collisions of other particles with 
the molecule, and different collisions can have dif- 
ferent effects. We shall assume that the probability 
distribution for these collisions in a particular in- 
terval, both for the times at which they occur and 
the effects which they have on the molecule, are sta- 
tistically independent of what happens in other in- 
tervals. Obviously this cannot be exactly correct, 
but on physical grounds it seems reasonable pro- 
vided At m is not too short, which is why we assume 
that it is larger than t c . In addition we assume, as 
is appropriate for a steady state situation, that this 
probability distribution depends only on the length 
At m of the interval, and not otherwise on t m . 

The next assumption is that at the beginning of 
a time interval of length At the molecule and the 
environment can be adequately described, for the 
purposes of what happens next, as a tensor product 
of a molecule state and some density operator for 
the environment. 1 The latter is the quantum analog 
of a probability distribution for incoming particles 
which might collide with the molecule during this 
time interval. This density operator can be "puri- 



Working out the general connection of decoherence with 
thermodynamic irreversibility and the properties of steady 
state is, at a fundamental level, an unsolved problem. Our 
hope is that its eventual resolution will justify present prac- 
tice by the experts in both the quantum information and 
decoherence communities, whose example we are follow- 
ing here. We note that some detailed quantum mechanical 
treatments of decoherence for situations similar to ours can 
be found in, e.g., [ToL 111! and references therein. 



fied" by regarding it as arising from an entangled 
pure state between the environment and an auxil- 
iary reference system, which we also take to be part 
of the environment. The overall time development of 
the molecule and the environment during the inter- 
val At is then given by a unitary time development 
operator, corresponding to an appropriate Hamilto- 
nian, acting on the system and environment, result- 
ing in an isomctry mapping the Hilbcrt space 7-Lm 
of the molecule onto Hm ®'He 1 where He is the 
Hilbert space of the environment. If one traces out 
the environment the result is a quantum operation 
or channel from H,m to itself: the channel input is 
the molecule at the beginning of the time interval 
At and its output is the molecule at the end of this 
interval. It is represented by a completely-positive 
trace-preserving (CPTP) supcroperator T(At) from 
the space Hm of linear operators on T-Lm to itself. 
Tracing out the molecule instead of environment at 
the end of the time interval results in a correspond- 
ing CPTP map T c {At) from Hm to He, represent- 
ing the complementary channel. See, for example, 
|12| for further details on how the direct and com- 
plementary channel are related to the isomctry. 

The assumption of statistical independence of suc- 
cessive time intervals, and that the environment is 
in a steady state, allows us to treat the interval from 
t\ to ti in the same way as the interval from to to t\ . 
Thus a succession of time intervals can be thought 
of, so far as the molecule is concerned, as a set of 
channels in series, with dynamics corresponding to 
an appropriate composition of the superoperators 
T{At m ) for the corresponding intervals. 

The use of a supcroperator T(At) that depends 
only on the length At of the time interval may give 
rise to the misleading impression that we are assum- 
ing exactly the same number and type of collisions 
for any interval of length At. But this is not so. To 
understand why, consider a classical stochastic pro- 
cess for which the independence of successive time 
intervals justifies using a Markov model, and for sim- 
plicity assume that all the time intervals are of equal 
length. Let j be a discrete index labeling molecule 
states at a single time, j m its value at the begin- 
ning of the m'th time interval, and M^ n \j',j) the 
Markov matrix for a transition j — > j' if precisely n 
collisions occur in one time interval. The probabil- 
ity distribution for a collection of histories that all 
begin in the state j\ conditioned on a specified set 
m, H2, • • -Tif of numbers of collisions in the different 
intervals is 

Pr(ji,j2,---j/+i 1 711,712,. ..n/) = 
M^(j f+1 ,j f ) ■ ■ ■ M^(j 3 , j 2 )M^\j 2 , n ) (3) 

On the other hand, if the sequence of collision num- 
bers is not known, the probability not conditioned 
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on this information is given by 

M {av) (j f+ uJ f ) ■ ■■M^(j 3 ,j 2 )M^(j 2 ,j 1 ) (4) 



where AT-/"- is the averaged Markov matrix, 
M^{j',j) = J2^(n)M^(j',j), 



(5) 



and Pr(n) the probability of n collisions during a 
single time interval. In the quantum case the single 
superoperator T(Ai) is the analog of an averaged 
Markov matrix, and it will allow us to correctly com- 
pute the probability of a sequence of histories as long 
as we do not try and condition it on more detailed 
information about the initial state of the environ- 
ment at the beginning of the time interval. 



B. Explicit form for the superoperator T 

As noted above, T(At) only makes physical sense 
for At greater than some correlation time r c . Keep- 
ing this in mind, it is nonetheless very convenient 
to think of the argument of T(At) as a continuous 
variable, which we shall hereafter denote by i, thus 
T{t). This superoperator can be written in various 
ways, e.g., using Kraus operators or as a matrix us- 
ing some basis of the operator space of a qubit. A 
convenient basis is provided by the Pauli operators: 
Co = I, <ti = X , oi = Y , 1T3 = Z, in terms of which 
we write 



(6) 



using a matrix T of real coefficients whose first row 
(because T is trace preserving) is (1, 0, 0, 0). The re- 
maining rows constitute a collection of 12 (real) pa- 
rameters which arc only constrained by inequalities 
that ensure that T is completely positive. Applying 
T(t) to a density operator p = Pj°~j at i = 0, 
with the coefficients {pj} regarded as a column vec- 



tor p, results in a density operator p = Pj&j at 
time t, where p = T • p. 

Rather than explore the entire parameter space, 
we have assumed that T has the particularly simple 
form 



T(t) 
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(7) 



where u) is the precession frequency for \R) to \L) 
and back again for the isolated molecule — the energy 
difference between |0) and |1) is Huj — and 7 > is the 
rate of decoherence. Justification based on scattering 



theory for this form of T has been discussed in 0, 

To see the motivation behind ([7| , first consider the 
case of the isolated molecule with no decoherence, 
7 = 0. Then 



T(t) = R(t) 



/l 0\ 

cosut — sinojt 

sinwi cos cot 

\0 1; 



(8) 



corresponds to precession about the z axis in a Bloch 
sphere picture. Next, suppose that uj = 0, so that 
only decoherence is present. Then 



T(i) « D(t) 
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(9) 



when t is small, and 



T{t) = D(i) • R(i) + Q(t 2 ) = R(f) • D(t) + 0(t 2 ). 

(10) 

with 0(t 2 ) a second order correction. Thus ([7]) com- 
bines the competing effects of decoherence and the 
molecule's internal dynamics. 

The motivation behind (J9J) is a simple physical 
picture in which if the environment is initially in the 
state \E) its interaction with the molecule during a 
collision corresponds to the unitary transformation 



\L) g> \E) -> \L) g> y\-2p\E) + 
\R) ® \E) -> \R) ® (v/l-2p|£:) + 



2p|^». 



(11) 



where the environment states and \E") arc 

orthonormal. The intuitive idea is that the distinc- 
tion between \L) and \R) is carried off to the distinct 
environmental states \E') and \E") with an ampli- 
tude that increases with p, a quantity lying between 
and 1/2 which is a measure of the effectiveness of 
the decoherence. The decohering effect is unchanged 
if \E) on the right side of (jTTJ) is replaced with any 
other state \E) as long as it is orthogonal to \E') and 
\E"), i.e., if the alteration does not depend upon the 
difference between \L) and \R). One can represent 
the channel corresponding to (fTl"]) by three Kraus 
operators associated with \E), \E') and \E"), but 
an equally good form uses just two Kraus operators 
^/l — p I and y/p X corresponding to a "bit flip" 
channel in [l~3[ p. 376. When p = there is no deco- 
herence (a perfect channel) whereas for p = 1/2 the 
collision "collapses" the molecule into either \R) or 
\L). If one setsp = the superoperator correspond- 
ing to the process pT) is given by © , and this makes 
sense for t of the order of the time between collisions. 
However, as noted above in Sec. Ill Al the superop- 
erator T{t) can appropriately represent a situation 
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in which the number of collisions in the interval t is 
a random quantity. 

The matrix T(i) and the density operator p(t) 
thought of as a column vector satisfy the simple lin- 
ear differential equations: 

The second is equivalent to a master equation in 
Lindblad form 

^ = -i[H >P }/h + 1 (X P X-p) > (13) 

with H/h = uiZ/2 as in (JTJ. 

Solutions to (fT2j) and (fT3|) can of course be ex- 
pressed as linear combinations of exponentials of the 
form e Xjt , where 

Ai = 0, A 2 = -7 + f, A 3 = A 4 = -2 7 (14) 



a = Vi 2 - ^ (is) 

are the eigenvalues of the matrix S. Note that A 2 and 
A3 occur in solutions of the form e A * to the damped 
oscillator equation d 2 x/dt 2 + 2 r ydx/dt + uj 2 x = 0. 
Thus for 7 < u they are complex conjugates of each 
other lying on a circle of radius u> in the complex 
plane, corresponding to oscillatory solutions, while 
for 7 > oj both are real and negative, correspond- 
ing to damped motion without oscillation. Critical 
damping 7 — lo corresponds to a phase transition in 
the sense of a changeover between two qualitatively 
different types of behavior. The explicit form of T(t) 
is given in Appendix |XJ 

III. GENERAL ASPECTS OF CONSISTENT 
HISTORIES AND INFORMATION 

A. Introduction to histories 

In the (consistent or decoherent) histories formal- 
ism a history is a sequence of quantum properties, 
identified by projectors onto appropriate subspaces 
of the quantum Hilbcrt space, at a succession of 
times ti < t 2 < • • • < */; see Ch. 8 of Q. In 
the situation at hand we use a sample space of mu- 
tually exclusive histories formed by assuming that at 
time t m the properties of interest to us correspond 
to a collection {P„ m } of projectors which form a 
decomposition of the identity: 

£P«'"=/, (P% r )i=P% r = (P% r f. (16) 

Here the subscript m labels the time, while the su- 
perscript a m is not an exponent but instead a label 



to differentiate the projectors at this time. Choos- 
ing at each time a property from the corresponding 
decomposition of the identity yields a history repre- 
sented by a projector 

Y a = P^ P 2 Q2 • • • Pj*, a = (ai,a 2 , ... a/) 

(17) 

on the history Hilbert space W = H ? = W0W0 
• • • H formed by the tensor product of the Hilbert 
space with itself / times. Here is a tensor prod- 
uct symbol with the same significance as 0, but em- 
ployed to distinguish different times. The physical 
significance of Y a can be seen by reading p7|) as 
"property P" 1 at time t\ followed by property P 2 ° 2 
at time i 2 followed by. . . ." 

For a closed system in which the unitary 
(Schrodingcr) time development from t m to t m +i 
is described by the operator U m +i tm , probabilities 
(probabilistic weights) can be assigned using the de- 
coherence functional [14j | 

B(Y a , Y?) = Ti[P" f U fJ -i ■ ■ ■ P.2 2 U 2 .iPi lx fo 

P i 1 U 1 , 2 Pi°---U f - 1J P^], (18) 

where ^0 is some initial state, provided the consis- 
tency conditions 

B{Y a , Y p ) = whenever a. ^ [3 (19) 

are satisfied. Here a 7^ (3 means that for at least 
one time t rn it is the case that a m 7^ j3 m . When 
(fT9| holds one assigns the positive weight W(a) = 
1D>(Y a ,Y a ) to the history Y a . The probability of 
each history is its weight divided by the total weight 
of all the histories; if ^0 is a normalized density 
operator this total weight is 1 and the probability of 
history a is W(a). 

The Hilbert space for the present discussion is 
H = Hm^He, where Hm is the Hilbert space of the 
molecule and He that of the environment. However, 
the histories of interest to us refer to properties of 
the molecule, not the environment, and we employ 
the usual convention that P"™ representing one of 
these properties can denote both a projector on Hm 
or its counterpart P% m Ie on H. For the initial 
state we let * = I M | <& E )(<& E |, where 

\* E ) = \Ex)®\Ez)®->\Ef-x) (20) 

is a "giant" tensor product state on the environment 
chosen in such a way that during the time interval 
between t m and t m+ i the molecule will interact only 
with the piece \E m ) in this tensor product in a man- 
ner determined by C/ m +i. m ; after that this part of the 
environment can be ignored so far as the molecule is 
concerned. In particular, if we take a partial trace 
over the environment of the middle portion on the 
right side of (fT5|) at time t 2 , the interaction of the 
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molecule with |£M is chosen so that 



Tr, 



U 2A P^oP^U 1 ,2 =T 2 .i(Pi 1 Pi 1 ), (21) 



where 7i, 1 = T(^2 — ^i) is the superoperator that 
maps the state of the molecule at the beginning of 
this time interval to its state at the end. In the 
same way, if the partial trace over the environment 
is carried out at time £3 the result will be 



TsAP^iiP^phPi 



(22) 



with 7js,2 = T(ts — £2), and similarly for later times. 
Consequently, for our model the decoherence func- 
tional is given by 



(y« Y P) = Tr M [P* f T f ,f-i(- ■ ■ P^TzAP? 



phPi 



(23) 



an expression which no longer makes any (direct) 
reference to the environment. See Sec. Ill of [l5| for 
a more detailed argument. 

If all the projectors in the decomposition {P r " m } 
are rank 1, which is to say they project onto pure 
states of the molecule, and the consistency con- 
ditions are satisfied, then the probabilities (corre- 
sponding to the diagonal elements of the decoher- 
ence functional (|23|) ) are those of a memory less hop- 
ping process - a Markov process. If the time steps 
t m are identical and the same decomposition is used 
at every time, this process is stationary (homoge- 
neous, i.e. same Markov matrix at each timestcp), 
but in general it is nonstationary (inhomogencous). 
Both cases are of interest for our model, as discussed 
below in Sec. IIVI 



B. Forwards and backwards conditions 

Finding collections of histories such that the con- 
sistency condition (|19| is satisfied is made somewhat 
easier by the following observation. Suppose it is 
the case that for every m between 1 and / — 1, if 
Q is a linear combination of the projectors in the 
set {Pm m }i then T m +i,m(Q) is a linear combination 
of the projectors in the set {Pm+V}- When this 
forward condition is satisfied, the family of histo- 
ries will be consistent, as can be seen in the fol- 
lowing way. The functional B(Y a , Y@) in ([23]) will 
vanish if a x ^ ft, since P^Pf 1 = 0. If the for- 
ward condition is satisfied, 72,i(Pi ai ) will be a linear 
combination of projectors in the collection {P^ 2 }, 
and will therefore commute with any projector in 
this collection. Consequently, P? 2 T2,i( p i 1 ) P 2 2 = 
P 2 Q2 P 2 & 72,i (Pf 1 ) will vanish whenever a 2 ^ ft 
since P^ 2 P^ 2 = 0, and if it docs not vanish it 



will be some linear combination of the {P^ 2 }- Pro- 
ceeding in the same way for larger to one sees that 
3(Y a ,YP) will vanish if, for any to, a m ^ /3 m . 

Note that the forward condition is a sufficient but 
not a necessary condition for consistency. The same 
is true of the backward condition: for every to be- 
tween / and 2 it is the case that if Q is a lin- 
ear combination of the projectors in {P^" 1 }, then 
7~m m -i(Q) is a linear combination of the projec- 
tors in {P,"™! 1 }- Here T^ im _ 1 denotes the ad- 
joint of the superoperator with respect to the Frobe- 
mus inner product: {V{A),B) = (A,T(B)) where 
(A,B) = TriA^B). The proof of consistency when 
the backward condition is satisfied proceeds in the 
same way as for the forwards condition, but in re- 
verse. Start with (|23|) and rewrite the argument in- 
side the trace by first cycling Py to become the first 
term, and then replacing T)j-i with Tj j_ 1 acting 
on Pj 1 Pj f , and continue this cycling process to con- 
vert all T to . [Note that since T is a (completely) 
positive superoperator, (J^{A))^ = T^{A^)] In the 
case of qubits, the situation of primary interest for 
the present paper, one can show that consistent fam- 
ilies of histories of the type (fT7)) must satisfy cither 
the forward or the backward condition. 



C. Measuring information 

We will want to discuss and quantify the informa- 
tion about the initial state of the molecule as time 
goes on. We can do this within the context of the 
histories formalism. Alternatively, we can do this in 
the context of the quantum channel formalism, i.e., 
quantifying the distinguishability of density opera- 
tors at the output of a quantum channel, and as we 
will see there is some connection between the two 
approaches. 

Let us first consider information from the histories 
perspective. Suppose some consistent family of his- 
tories uses the projective decompositions Pi = {P(} 
at time t\ and P m = {P„} at time t m , with t\ < t m 
[for simplicity here we replaced the indices ct\ and 
a m in (fT6]) and (fT7]) with j and k]. As in fl6j ]. we 
will equate the notion of a projective decomposition, 
like Pi , with a type of information about the system, 
in our case the molecule. A convenient measure of 
how much of the Pi type of information about the 
molecule remains at time t m is the Shannon mutual 
information 



P(P :P m ) = H{P x )+H(P m ) - H(P u P m ), (24) 

where H(P\) is the familiar Shannon entropy. In 
particular if Pi, P m , and P m < are projective decom- 
positions associated with a consistent family at three 
successive times, and if the probabilities correspond 



7 



to a Markov process, then (see, e.g., p. 510 of (l3j ) 
if (Pi : P m >) cannot be greater than if (Pi : P m ): the 
information about the initial situation can only de- 
crease with time. For simplicity, in what follows 
we will set Pi(Pf) = 1/di for all j, where d\ is the 
number of projectors in the decomposition Pi. Then 
H(Pi : Pi) = H(Pi) = logdx, and hence the infor- 
mation decays from its initial value of log di as time 
goes on. 

Now, alternatively, consider the quantum chan- 
nel perspective, where we will quantify how much 
of the Pi type of information remains at time t m 
by measuring the distinguishability of the condi- 
tional density operators at the output of the rele- 
va nt q uantum channel. (This approach was taken 
in [12j.) To measure distinguishability of density 
operators, in particular if these density operators 
do not commute, we need a measure that is inher- 
ently quantum-mechanical, which is provided by the 
Holevo function 

x({Pj,Pj}) ■= S(Y^PjPj) - ^ZP]S{pj) (25) 

3 3 

defined for an ensemble {pj , Pj}, where pj is the 
probability assigned to the density operator pj, 
and S(p) := — Tr(plogp) is the von Neumann en- 
tropy. Applying this measure to the ensemble 
{1/di, T m .i(P()}, where T m ,i is the quantum chan- 
nel that governs the molecule's evolution from ti to 
t m , gives a quantitative measure of how much Pi 
information remains at time t m , and we write this 
as 

x(p l! r W ) = x(p l! r m , 1 ):=x({^^^}), 

(26) 

where t = t m — t\. 

Equations ([24]) and (|26|) give two alternative ways 
to measure the loss of information from the sys- 
tem over time. Equation (|24|) has the advantage 
of a clear conceptual interpretation, whereas Equa- 
tion (|26p has the advantage of being easy to compute 
since one does not need to go through the histories 
analysis to compute it. Fortunately, there is a con- 
nection between these approaches. It turns out, see 
the argument in Appendix [D] that for a family sat- 
isfying the forward consistency condition 

H(P 1 :P m )=x(Pi,T m ,i). (27) 

A similar sort of connection holds for families sat- 
isfying the backward consistency condition (but in- 
volving the adjoint channel 7^ i), but for simplicity 
we will focus on families satisfying the forward con- 
dition to illustrate information flows in Sect. [V] 

One can also quantify information flow from the 
molecule to the environment with the quantum 



channel approach by using the complementary chan- 
nel with superoperator T c , introduced in Sec. Ill Al 
In fact, T c is completely determined by T up to 
an isometry on its output (the environment), which 
does not affect distinguishability measures like x- 
Hence, the following information measure is well- 
defined: 

X(Pl,T C (t)) = X (Pl,Vn A ) := X ({ ^, ^7^»: 

«i Ir(Pj') 

(28) 

where t = t m — t±. It quantifies the amount of 
the Pi type of information about the molecule (at 
time ti) that is present in the environment at time 
t m . Though one cannot in general equate this 
x(Pi,7^j i) with a Shannon mutual information be- 
tween the molecule and the environment, the former 
provides, as is well-known (e.g., p. 531 of [IH), an 
upper bound on the latter. 

We note that there can be a tradeoff in sending 
information to the environment and preserving it in 
the molecule, which is most dramatic for comple- 
mentary or mutually-unbiased bases Pi and P{: 

x(Pi,T m ,i) + x(P[,T^i) < logdi. (29) 

This inequality is from Corollary 6 of [l2j ]. 

IV. CONSISTENT FAMILIES FOR OUR 
MODEL 

A. Differential equations 

Our model has only two states, and therefore any 
(nontrivial) decomposition of the identity involves 
only projectors of rank 1 onto pure states. Thus 
a consistent history family corresponds to a two- 
state Markov process (sometimes called a "telegraph 
process"); in general this process is nonstationary: 
the transition rates depend upon the time. While 
such a process can be discussed using discrete times 
separated by finite intervals, the results are simpler 
and the mathematical expressions more transparent 
if one adopts a continuous time approximation with 
differential equations in place of difference equations. 
It should, of course, be kept in mind that the pro- 
cesses here described are not truly continuous, since 
time intervals shorter than the correlation time r c 
introduced in Sec. HTl lack physical significance. The 
continuous time approach should be satisfactory as 
long as both lot c and 7T C are small compared to 1. 
Note that this condition can still be true even when 
7 is large, as long as r c decreases as I/7, which seems 
physically plausible. 

For families satisfying the forward consistency 
condition the relevant differential equations can be 
obtained in the following way. At a particular time 
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the decomposition of the identity will correspond to 
two projectors, call them po and p\, represented by 
end points or antipodes of a diameter of the Bloch 
sphere. Let the direction of this diameter be de- 
noted by the usual polar and azimuthal angles 9 and 
0: the z axis at 6 = and the x axis at 9 = 7r/2, 
= 0. Which end of the diameter corresponds to 
these angles does not matter for the following discus- 
sion. The locations of these end points after a short 
time interval is determined by the master equation 
• One can show that because of the form of S in 
([7]) they are still located on a diameter of the Bloch 
sphere, but are now a bit closer to its center. The 
rate of change of the diameter's direction is repre- 
sented by the differential equations 

d(f> d9 o 

— = lo — 7 sin 20, — =7 sin 20 cos 0, (30) 

whereas the shift towards the center can be used to 
calculate the instantaneous transition rate 

K = 7(l-sin 2 0cos 2 0), (31) 

which enters the rate equations 

dpo/dt = k(-pq +pi), dpi/dt = K(p -pi) (32) 

for the probabilities associated with these two states. 
The backwards consistency condition can be ana- 
lyzed in a similar way, and leads to the differential 
equations 

dd> d9 2 

— = cj + 7 sin 20, — = — 7 sin 20 cos 0. (33) 

dt dt 

governing the direction of the diameter, and to ex- 
actly the same expression pip for the transition rate. 
For a more detailed derivation of these formulas see 
Appendix [B] 

B. Stationary families 

If the angles 9 and which determine the diam- 
eter for the projectors forming a consistent family 
do not change with time the Markov process is sta- 
tionary or homogeneous, in the sense that the states 
and the transition probabilities do not change with 
time; of course the actual state of the molecule is 
varying randomly as it hops back and forth between 
the two states. The simplest case is what we call the 
z family, in which 9 = (or 7r), thus d9/dt = in 
(|30f or ([33| and d(f>/dt is irrelevant. The two projec- 
tors (I + a z )/2 and (I — a z )/2 correspond, respec- 
tively to the even parity (ground) and odd parity 
(excited) states of the isolated molecule. Thus we 
have a two-state stationary Markov process in which 
the molecule spends a certain amount of time in the 
even parity state before flipping instantaneously (in 



our continuous time approximation) to the odd par- 
ity state where it remains for a random time inter- 
val before flipping back. The time r between flips is 
a random variable with an exponential distribution 
e~ 7T , since setting 9 = in (j3Tj) gives 

« 2 = 7 (34) 

for the transition rate. On average the molecule 
spends an equal amount of time in both states, which 
means that in our model the environment has an ef- 
fective temperature T ^> hio/kB- 




FIG. 1: (a) The steady-state solutions for correspond 
to the intersections of sin 20 (solid curve) with 
(dashed lines), shown here for 7 > uj > ("strong deco- 
herence" regime), (b) These steady-state solutions are 
plotted schematically on the Bloch sphere, as if the z- 
axis is going into the page. 

In addition to the z family just discussed there 
are stationary families in which the projectors cor- 
respond to points in the x-y or equatorial plane of 
the Bloch sphere, so 6 = -k/2 with d(f>/dt = in ([30} 
and (|33|) . and thus 

sin 20 = ±(w/7). (35) 

For < cj/7 < 1 there are four solutions as shown 
in Fig. [TJ which coalesce into two for lu/j = 1. For 
lj/j > 1 these families disappear, leaving the z fam- 
ily as the only stationary family. In the limit of 
strong decoherence, small w/7, two of the families 
approach the x axis and two the y axis of the Bloch 
sphere, so we shall refer to them as the dressed x- 
and dressed y-familics. The associated transition 
rates k x and n y are given by —1/2 times the cor- 
responding eigenvalues of S, see (fl4]) : 

Kx = -A a /2 = (7-0/2, Ky = -A3/2 = (7 + O/2. 

(36) 

As 7/w becomes very large the dressed a;- families 
approach the x or chirality basis \R) and \L) of @, 
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and the transition rate k x « lj 2 /At/ becomes very 
slow. Thus these families represent long-lived (al- 
most) chiral states when decoherence is rapid com- 
pared with the the tunneling rate. (But see the fur- 
ther discussion in Sec. II V EH 



C. Nonstationary families 

The equations (|30|) and (f3"3")l can be integrated 
in closed form to obtain the bases corresponding 
to nonstationary consistent families, Appendix [C] 
However, the solutions are fairly complicated expres- 
sions. The time evolution for some cases in which 
lj = 1 and 7 is either less than or larger than 1 is 
shown in Fig. [2l For 7 < u) the diameter rotates 
continuously about the z axis (the discontinuities in 
cj) are of course artifacts of the plot) with an angular 
frequency of 



V 



(37) 



while the polar angle 9 tends either to ir/2 for the 
forward or to (cquivalcntly, ir) for the backward 
consistency condition. In the limit in which 7 tends 
to 0, no decoherence, one has a simple rotation of 
the diameter of the consistent family about the z 
axis at a rate lj with 6 fixed. This same tendency 
is seen in the dependence of 9 on time when 7 > lj, 
whereas <j> more or less rapidly approaches one of the 
values corresponding to a stationary family. Note 
that along with the continuous change of basis there 
is a random flipping from one of the basis states to 
the other at a rate given by (|3"Tj) . so one is dealing 
with a nonstationary Markov process. 

Thus in the Bloch sphere picture, for 7 < lj, 
the families "corkscrew" about the z axis (going 
away from or towards this axis for the forward or 
backward families, respectively), with 4> periodically 
coming back to the same value at time intervals that 
are integer multiples of tt/t). If 6 = 7r/2 it remains 
constant, so the same basis reoccurs after an inter- 
val of 77/77. ^ om y these discrete times are consid- 
ered, the result is what one might call a stroboscopic 
family which can be thought of (Xti (X discrete time 
stationary Markov process. 



D. Phase transition 

As the parameters 7 and lj vary there is a phase 
transition, a qualitative change of behavior, when 
they are equal. This manifests itself in a variety of 
related ways. For 7 < uj the eigenvalues of S include 
a complex-conjugate pair A2 and A3, (|14[) . which coa- 
lesce into a single degenerate eigenvalue at the tran- 
sition, and thereafter, for 7 > lj, become a pair of 
distinct real eigenvalues. This is, of course, precisely 



the behavior one finds in a classical one-dimensional 
oscillator when the damping passes through the crit- 
ical value. For 7 > lj these eigenvalues are the de- 
cay rates for the dressed- a; and dressed-y continuous 
stationary Markov processes discussed in Sec. IIVBI 
On the other hand, as 7 decreases towards lj from 
above, the four stationary families shown in Fig.fljb) 
coalesce into two, corresponding to diameters of the 
Bloch sphere midway between the x and y axes, and 
for 7 < lj they no longer exist: the only remain- 
ing continuous stationary family is the z family. As 
noted above in Sec. IIV CI for 7 < lj there is a new 
class of "stroboscopic" families defined using a pe- 
riodic time interval. As 7 approaches lj from below 
this period becomes infinitely long. The behavior 
of nonstationary continuous families is also different 
for 7 < lj and 7 > lj. For the former (j) increases 
indefinitely and monotonically with time, although 
this motion, which is simply linear when 7 = 0, be- 
comes more and more "jerky" as 7 increases towards 
lj. See the example for j/lj =1/2 in Fig. [2j For 
7 > lj, <j) approaches a fixed value with increasing 
time, and no longer "winds." (Again, the damped 
harmonic oscillator provides a helpful analogy.) One 
might say that the nonstationary continuous families 
transition from a damped oscillatory character to a 
purely damped character as 7/w increases, passing 
through the critical value of 1. 

In terms of its mathematical structure as repre- 
sented in the master equation this is a dynamical 
quantum phase transition of the sort discussed in 
quantum optics for two level systems in Ch. 11 of 
17 1 and in [l8|, and in a more general context in 
It appears that the vanishing of the inver- 
sion transition in ammonia is of this type; see the 
discussion in Sec. IVIBl below. We believe that ours 
is the first attempt to explore dynamical properties 
near such a transition using the histories approach. 



E. Physical interpretation 

Each consistent family contains a collection of his- 
tories, and each history a particular succession of 
micrscopic properties (subspaces of Hm)- One and 
only one history from this collection will describe 
the behavior of a particular molecule during a par- 
ticular interval of time. There is no need to make 
any reference to measurements, though it is in prin- 
ciple possible (i.e., does not violate the laws of quan- 
tum mechanics) to use a succession of suitably ideal- 
ized measurements to determine which of these his- 
tories is actually realized. But because one is dealing 
with a system exhibiting "quantum" behavior, i.e., 
in a regime in which a classical description is not 
adequate, it is important to keep in mind certain 
respects in which quantum descriptions differ from 
their classical counterparts. 
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FIG. 2: Time evolution of the Bloch-sphere angles <f> and 8 for the consistent description associated with the (a) 
forward and (b) backward conditions. In this case, </>(0) = 0, 8(0) = 0.2, uj = 1, and 7 = 0.5, 1.2,4 respectively for 
the solid, dashed, and dotted curves. The angles <j> an< l # are m units of radians, and the time t is in units of l/co. 



In particular, two consistent families of histories 
may be mutually incompatible with each other in 
such a way that they cannot be combined into a 
single description that makes sense. A well-known 
example is a spin-half particle where incompatibility 
arises from the fact that the operators for angular 
momentum in different directions do not commute 
with each other, and hence have no common eigen- 
vectors. It makes good (quantum) sense to say, for 
example, that S x = +1/2 (in units of ft), or that 
S z = +1/2, but there is no quantum property, no 
subspace in the Hilbcrt space, that corresponds to 
S x = +1/2 and S z = +1/2. So one cannot as- 
cribe simultaneous existence to S z and S x . In the 
consistent histories approach this inability to com- 
bine incompatible descriptions is codified as the sin- 
gle framework rule, and the consistency conditions 
discussed above in Sec. IIIII serve to extend this rule 
from a single time to a sequence of times. In ad- 
dition, just as two incompatible consistent families 
or descriptions cannot be combined into a single de- 
scription, they also cannot be compared: it makes 
no sense to ask which of two incompatible families is 
the "correct" one, or to look for some law of nature 
that single out one against another. Each consis- 
tent family provides its own quantum description in 
a way roughly analogous to looking at a mountain 



from different locations. For a detailed discussion of 
these points we refer the reader to [22T |. 

With reference to a tunneling molecule, consider 
the situation in which 7 is much larger than w, 
strong dccohcrence. There is a stationary family, 
Sec lIVB] the "crural" family, in which the molecule 
hops back and forth at a comparatively slow rate 
between the (dressed) left-handed and right-handed 
chiral states. (There are actually two of these 
dressed- x families, but when decoherence is strong 
there is very little difference between them.) This is 
the family to use if one is interested in understanding 
why a specific chirality, the left or right-handed form 
of the molecule, can persist for a very long time in 
a situation of strong decoherence. It provides a de- 
scription in terms of a stochastic two-state Markov 
process in which the rate of hopping from left to 
right-handed or vice versa is a well-defined function, 
k x in (|36p , of the parameters that enter the model. 
(Each hop is instantaneous on the time scale used for 
our description, in which intervals less than the cor- 
relation time t c do not enter; see Sec. [ill) As the rate 
of decoherence decreases, the hopping time becomes 
shorter and the amount of "dressing" required to 
produce a consistent family increases, which means 
that even though this family continues to provide 
a correct quantum description, it no longer corre- 
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sponds to a simple physical picture of a definite left- 
or right-handed molecule when 7 becomes compara- 
ble to OJ. 

In addition to this chiral family there is a z or 
"parity" family in which the molecule hops back and 
forth at random time intervals between parity eigen- 
states (energy eigenstates of the isolated molecule) , 
at a rate given by the decoherence rate 7, see (|34j) . 
The parity family is incompatible with the chiral 
family discussed above, and they cannot be com- 
bined. One should not try and imagine them as 
going on simultaneously; to do so would be to make 
the same mistake as supposing that S x and S z for a 
spin-half particle can simultaneously possess values. 
On the other hand, just as it is possible to mea- 
sure either S x or S z , but not both simultaneously, 
it is also possible in principle (without violating the 
laws of quantum mechanics) to determine by mea- 
surements the succession of events that occur in a 
parity family, or by a different set of measurements 
those occurring in a chiral family. Thus a relatively 
rapid but random flipping back and forth between 
parity eigenstates is a valid physical picture of the 
succession of microscopic states of the molecule, one 
which can be used both when the decoherence is 
strong and when it is weak. There is in addition a 
third stationary family for 7 > u, the dressed y fam- 
ily, which has a relatively rapid hopping rate in the 
strong decoherence regime, and eventually merges 
with the chiral family as the decoherence rate de- 
creases. We do not have a simple name or physical 
interpretation for this family. 

In the regime where decoherence is weak, 7 < oj, 
there are no truly stationary families, apart from 
the parity family discussed above. A relatively sim- 
ple nonstationary family is the one that employs 
an "equatorial" basis in the x-y plane of the Block 
sphere, 6 = tt/2, rotating at an average angular 
speed 77, see (|37|) . Let us call this the "tunneling" 
family, since it corresponds in physical terms to the 
molecule oscillating back and forth between the two 
potential wells. As 7 increases the rate of tunneling 
decreases and eventually goes to zero at the phase 
transition 7 = uj. In addition to the tunneling, the 
phase 4> undergoes random changes by 7r, instanta- 
neous on the time scale we are using, at a rate, (|31l) . 
proportional to 7, but also depending on the value 
of (f>. Thus we have a nonstationary Markov pro- 
cess. The random flipping rate increases with 7 at 
the same time as the tunneling rate is decreasing, so 
the simple physical picture of the molecule tunnel- 
ing from one potential well to the other breaks down 
upon approaching the phase transition 7 = uj. For 
larger values of 7 this consistent family no longer 
exists. 



V. INFORMATION FLOWS 

In the previous section we found various consis- 
tent frameworks for discussing the stochastic trajec- 
tory (for our model) of a tunneling molecule. We 
now wish to study the dynamics of information, e.g., 
the loss of information about the molecule's original 
state as time progresses. Section lTlI Cl discussed how 
the Shannon mutual information between the orig- 
inal state and the state at some later time, for the 
forward or backward consistent family, is equivalent 
to a particular Holevo \ quantity. Here, for simplic- 
ity, we will focus on families satisfying the forward 
condition, for which the information remaining in 
the molecule is given by (|27|) . and that flowing to 
the environment by (|28p . 

These quantities are shown in Fig. [3] for P\ = Z 
(parity basis) and P\~X (chirality basis) , both for 
information remaining in the molecule T{t) and that 
flowing to the environment T c (t). Figure [S^a) shows 
a case of strong decoherence, -f/u> = 2.5, where the 
curves as a function of time are quite smooth, con- 
sistent with the fact, Fig. [Ha), that the consistent 
family is rapidly approaching a stationary family. 
For weak decoherence, Fig. [3f b) with j/oj = 0.05, 
the consistent family is not stationary and the alter- 
nating rises and plateaus reflect this fact. The top 
curves in both (a) and (b) represent the sums, see 
(|39jl below, for one type of information remaining in 
the molecule and a mutually unbiased type flowing 
to the environment. 

For short times the individual information mea- 
sures can be computed, using the expressions 
for T and T c in Appendix [X] and noting that 
5(T C (|^)(?M)) = S'(T(|V>(V'I)) for any pure state 
IV'), to obtain 

x(X,T(t)) = l-0[t 2 \og(l/t)}, 
X(Z, T{t)) = 1 - <yt log(l/ 7 i) - 7* + 0[t 2 log(lA)], 
x(X,T c (t)) = 7*Iog(l/7*) + It + 0[t 2 log(lA)], 
x(Z,T c (t))=0[t 2 log(l/t)], (38) 

Here 0[ ] means that the correction term is of this or 
possibly some higher order. These expressions are 
consistent with x(A, T{t)) and x(Z, T c {t)) having 
zero slope at t = 0, and x(Z, T{t)) and x(X,T c {t)) 
having infinite slope at t = 0, as depicted in Fig. [3] 
That x(Z,T c (t)) in (|38|) has no term linear in t 
seems plausible in that the decoherence mechanism 
in our model has been chosen specifically to carry X 
information into the environment. The uppermost 
curves in Figs.[3]Ja) and[3jb) represent 

x(Z,T(t)) + x (X,T c (t)) = 

x(z,r c (t)) + x (x, T{t)) = 1 - o[t 2 iog(iA)], 

(39) 
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FIG. 3: Information flows in terms of the \ information 
measure (see text for definition), (a) in the strong deco- 
herence regime with 7 = 2 and ui — 0.8, and (b) in the 
weak decoherence regime with 7=1 and uj = 20. The 
measure \ is in units of bits, and the time t is in units 
of 0.8/cj in (a) and 20/a; in (b). 



where the first equality comes from Theorem 3 of 
[lH, and the second from (|3"5)) ; the correction must 
be negative in view of the bound in (|29l) . 

Another consequence of (|38|) is that for small t 
the flow of chiral information to the environment 
is compensated by a decrease of parity information 
remaining in the molecule: 



: x(X,T c {t)) = --x(Z,T(t)) + 0[t\og(l/t)}. 

(40) 

However, as noted above, both sides of (|40|) diverge 
logarithmically as t — > 0. Of course, these expres- 
sions lack physical meaning for times shorter than r c , 
and thus the divergence is a mathematical artifact. 
Nonetheless, this makes it difficult to define rates of 
flow of information in a mathematically clean way 
using the x measure. 

An alternative which avoids the divergence is to 
replace the von Neumann entropy S in the definition 
(|25p with the quadratic entropy 



S Q (p) = l-Tt(p 2 ). 



(41) 



In the case of a qubit channel with W an orthonor- 
mal basis, projectors W 1 + W 2 = I, and Pauli oper- 



ator <j\y :— W 1 — W 2 , the measure defined in (f2"rj|) 
becomes 



XQ {W,T(t))= l -Tv 



(tm)' 



(42) 



with a similar expression for the complementary 
channel if T is replaced by T c . One can use the 
usual Pauli representation to write 



aw 



n ■ rr. 



(43) 



where the x, y, and z components of the aw Pauli 
operator are given by the unit vector 

n = {n x -,n y ,n z } = {sm6 cos </>, sin 9 sine/), cos#}. 

(44) 

Using the short-time expressions for T and T c given 
in Appendix [A] one finds that for the direct and 
complementary channels 

X Q (W, T{t)) = 1 - Ayt [1 - n 2 ] + 0(t 2 ), (45) 



XQ (W,T c (t))=A 1 tn 2 x + 0(t 2 ). 



(46) 



Setting W = X in ((46]) and W = V in (|45]), where 
V is some basis in the y-z plane and thus mutually 
unbiased relative to X, the analog of (|40|) with x 
replaced by x,Q is 



dt 



X Q (X,T c (t)) 



dt 



XQ(^T(t))=4 7 (47) 



at t = 0, so the derivatives are now finite. Thus 
when one uses the XQ measure the rate of flow of 
X information to the environment equals the rate of 
decrease within the molecule of any type of informa- 
tion associated with a basis in the y-z plane. 



VI. D 2 S 2 AND NH 3 
A. D 2 S 2 

An order-of-magnitudc calculation of the decoher- 
ence rate 7 for a D2S2 molecule immersed in a gas 
of helium atoms shows that this system is in the 
strong decoherence regime under typical conditions. 
The flux q of helium atoms (atoms per unit area per 
unit time) is their concentration times their average 
velocity, given by q = P ■ ^8/ (irmnckbT), assum- 
ing an ideal gas with pressure P and temperature 
T, with mHe and kb the mass of a helium atom and 
Boltzmann's constant [2j|. At room temperature 
T = 300 K and a pressure P of 1 atmosphere this 
gives q w 3 x 10 28 atoms s _1 m~ 2 . Multiplying q 
by approximate cross sections a co \ 



1000a 2 , and 



Cdec ~ 100a 2 , for collisions and decoherence, taken 
from Q, with ao the Bohr radius, leads to a collision 
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rate of approximately 9 x 10 10 s 1 and a decoherence 
rate of 

7K9X10V 1 . (48) 

The estimated tunneling rate is w « 176 rad/s after 
correcting 2 the value published in Q by a factor of 
2tt. Thus 

7/cj « 5 x 10 7 , (49) 

which means strong decoherence, for which our re- 
sults in Sec. IIVI indicate that chirality is both a 
consistent description as well as a stable property, 
to a very good approximation. Probing the regime 
7 < uj for this molecule would seem quite difficult 
as it would involve very low pressures. Replacing 
deuterium with hydrogen and/or sulfur with oxygen 
leads to chiral molecules, e.g. H2O2, that have signif- 
icantly higher tunneling frequencies p4| and hence 
may be candidates for probing the 7 < uj regime. 

B. NH 3 

In the electronic ground state of ammonia NH3 the 
nitrogen lies to one side of the plane defined by the 
three hydrogens, but there is a relatively low poten- 
tial barrier separating it from the mirror image state 
on the other side, and tunneling in this double-well 
potential exhibits itself in the well-known inversion 
transition at a frequency of about 24 GHz. Although 
the molecule is not chiral, when it is rotating about 
an axis passing through the nitrogen and the mid- 
point between the hydrogen atoms the symmetry op- 
eration of parity (inversion) moves the nitrogen to 
the other side of the plane, changing the sign of the 
electric dipole, while leaving the angular momentum 
unchanged. Consequently, the energy levels with a 
nonzero quantum number K for this component of 
angular momentum arc split into two parity eigen- 
states by the inversion transition in a way similar to 
that in a chiral molecule. 

The tunneling transition has been observed di- 
rectly by microwave absorption, which at low pres- 
sure exhibits a set of closely-spaced lines associated 
with the different rotational states [25[ . As the pres- 
sure increases the lines broaden and merge, and the 
center of the merged line shifts towards lower fre- 
quencies, reaching zero frequency at a pressure of 
about 2 atmospheres (2(|. It has been suggested, 
e -g- 113 > that at pressures above this transition the 
ammonia molecule adopts a "pyramidal" shape with 
the nitrogen on one side of the hydrogen plane, anal- 
ogous to the shape of a chiral molecule with a def- 
inite handedness. Deuterated ammonia ND3 shows 



2 Private communication from K. Hornberger 



similar behavior, except that the low pressure tun- 
neling frequency is now at 1.6 GHz, and the center 
of the broadened line tends to zero frequency at a 
pressure of 0.12 atmospheres J28|. 

The shift towards zero frequency has been ana- 
lyzed theoretically using two different approaches. 
The first, exemplified by (29| , and with similar ideas 
in [lO, HH among others, starts with the observa- 
tion that since the ammonia molecule possesses a 
significant electric dipole moment when the nitro- 
gen is on one side of the hydrogen plane there will 
be a strong dipolc-dipolc interaction between nearby 
molecules. It is then proposed that this produces a 
sort of mean-field effect in which the polarization 
of one molecule influences its neighbors in such a 
way that eventually as the pressure increases and 
the molecules come closer together, the double well 
potential for a single molecule is changed into one 
with a single minimum on one side of the hydrogen 
plane, resulting in molecules of pyramidal shape. 

An alternative approach found in [32I - I35I ] focuses 
instead on the decohering effects of collisions be- 
tween gas molecules. It is argued that these colli- 
sions in addition to broadening the lines can also 
lower the tunneling frequency as the pressure, and 
thus the collision rate, increases. From this per- 
spective the electric dipole-dipole interaction, while 
significant in determining the collision cross section 
and the effects of collisions, is not the fundamental 
source of the line shift to lower frequencies. The lat- 
ter ought still to be present if ammonia is a dilute 
component in a nonpolar buffer gas. Of particu- 
lar significance for this second point of view is the 
work of Ben-Reuven [32|, H|| , who argued on theo- 
retical grounds that when proper account is taken 
of the effects of collisions the line shape, absorption 
as a function of frequency, is not adequately repre- 
sented by the Van Vleck and Weisskopf formula [36[ 
used earlier in [26j to analyze the experimental data. 
He proposed an alternative line shape function with 
three parameters, 7, £ and 5, proportional to the 
collision rate, and thus the pressure, to fit the ex- 
perimental microwave absorption data for NH3 and 
(with a different choice of parameters) ND 3 over a 
range of frequencies and pressures sufficient to in- 
clude that at which the tunneling frequency goes to 
zero. It is noteworth that this fit was achieved for 
all pressures and frequencies using just these three 
parameters, whereas the earlier analysis of Bleaney 
and Loubser [26| was carried out by adjusting two 
parameters separately for each pressure. 

The validity of Ben-Reuven's analysis is supported 
by the fact that more recent data on microwave ab- 
sorption by ammonia in mixtures of hydrogen and 
helium (of interest in studies of the atmospheres of 
Jupiter and the other giant planets) has been fitted 
usin g hi s line shape formula for the tunneling transi- 
tion [37J with, of course, different choices of paramc- 
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ters for the different species scattering from the am- 
monia molecule. Since neither hydrogen nor helium 
has an electric dipole moment, this tends to support 
the idea that collisions, rather than dipole-dipole in- 
teractions as such, are what drive the transtion fre- 
quency to zero in pure ammon ia g as as the pressure 
rises. The numbers given in [371 ] would suggest a 
phase transition at about 20 atmospheres for ammo- 
nia in a buffer gas of hydrogen at room temperature. 
Replacing NH3 with ND3 should bring the transition 
pressure down by a factor of 15, and replacing hy- 
drogen with a gas of some other nonpolar molecule 
might be advantageous. Thus a direct experimental 
test of whether dipole-dipole interactions are or are 
not essential for understanding the vanishing of the 
tunneling frequency seems feasible. 

Our very simple decoherence model corresponds 
to setting 7 = ( and S = in Bcn-Rcuven's theory as 
it applies to a two-level system. In fact, he achieved 
a good fit to the experimental data with 6 = 0, but 
with 7 larger than £ by a factor of around 1.3, see 
p. 21 of [33I . To have 7 larger than £ in our model 
would require our adding another source of decoher- 
ence. The phase transition present in our model is 
also clearly present in Ben-Reuven's work; see the 
discussion of the spectrum of the perturbed Liou- 
ville matrix in Sec. 4C of [H|, where the eigenvalues 
change character when £ passes through the value 
ujq + 5; this is the counterpart of our 7 = 0;. Hence 
it seems that the vanishing of the inversion frequency 
in ammonia with increasing pressure is an instance 
of the sort of phase transition that occurs in our 
model. A more detailed comparison, which we have 
not attempted, would require our including an addi- 
tional mechanism for decoherence to make 7 larger 
than and dealing with complications caused by the 
presence in ammonia of a number of different rota- 
tional states. Nonetheless, we think our considera- 
tions provide some insight into the sense in which the 
ammonia molecule can be said to be "pyramidal" in 
the gas at high pressure and lack this feature at low 
pressures. Namely, when collisions are sufficiently 
frequent there is a consistent family of histories, the 
chiral or x family discussed in Sec. IIVEI in which 
in quantum mechanical terms the molecule spends a 
time much longer than the tunneling time in a pyra- 
midal shape (or a "dressed" state close to it) which is 
"chiral" in the sense that electric dipole moment has 
a definite orientation relative to the angular momen- 
tum, with occasional random hops between the two 
pyramidal possibilities. As the pressure decreases 
towards the transition pressure the pyramidal pic- 
ture begins to break down: the hops become more 
frequent between dressed states, which are starting 
to lose their pyramidal character. At still lower pres- 
sures it is better to think of the molecule as contin- 
uously tunneling back and forth, rather than pos- 
sessing a fixed pyramidal form, with a period that 



diverges as the pressure rises to its value at the tran- 
sition. 



VII. CONCLUSION 

We have shown how the decoherence of a two-state 
tunneling molecule, a chiral molecule or ammonia, in 
the presence of a buffer gas can be described in terms 
of a succession of quantum states of the molecule 
itself that form a consistent family of histories, on 
a sufficiently coarse time scale so that intervals arc 
always longer than a correlation time. Our model is 
described by just two parameters, a tunneling rate w 
and a decoherence rate 7, and its essential properties 
depend upon the ratio j/u>. In addition we have 
studied the flow of information to the environment, 
along with its retention by the molecule itself, during 
the process of decoherence. 

We found a large variety of consistent families, 
some of which are stationary (in the sense of Markov 
processes) and some of which are not. In the regime 
7/cj ^s> 1 of strong decoherence there is a station- 
ary family (actually two closely related families) in 
which the molecule spends a relatively long time in 
one of its chiral states before flipping to the other 
chirality, and eventually flipping back again, in a sta- 
tionary Markov ("telegraph") process, with a tran- 
sition rate which is approximately w 2 /47 for 7 ^> ui, 
and hence quite slow compared to the tunneling rate 
lj. Thus this "chiral" family explains the persistence 
of chirality for a long period of time when there is 
strong decoherence. However, as 7/w decreases, the 
transitions between chiral states become more fre- 
quent and the states themselves (the "dressed x" 
states of Sec. HVBj) lose their chiral character, un- 
til finally this family disappears entirely at a phase 
transition 7/w = 1. 

We have found two other stationary consistent 
families for j/us > 1. One of them is the "parity" 
family in which the molecule is at each of the times 
considered in one of the two states of definite par- 
ity (the energy eigenstates of the isolated molecule) , 
but with a transition rate of 7 between them, thus a 
rapid flipping compared to transitions between chi- 
ral states when the decoherence is large. This family, 
present at all values of 7/w, is incompatible, in the 
quantum mechanical sense, with the chiral family: 
while both provide valid quantum descriptions, they 
cannot be employed simultaneously; see the discus- 
sion in Sec. IIVEI The other stationary family for 
7/u; > 1 (again there are actually two families) is 
the "dressed y" family of Scc. UVBl It involves a rel- 
atively rapid flipping between two orthogonal quan- 
tum states for which we do not have a simple phys- 
ical interpretation. Like the chiral family this one 
only exists for 7 > uj. 

For 7 < uj there is a nonstationary "tunneling" 
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family in which the molecule oscillates back and 
forth between the two chiral states at a rate that 
goes to zero as j/u increases to 1. On top of this rel- 
atively smooth oscillation there are random changes 
in phase which constitute a nonstationary Markov 
process, with a rate that increases with 7. This 
tunneling family disappears at the phase transition 
7 = 0;. The only truly stationary family in the 
regime 7 < ui is the parity family. In addition, both 
for 7 < uj and for 7 > uj there are a variety of non- 
stationary consistent families in which the orthogo- 
nal basis used to describe the quantum system tends 
with time towards one of the stationary families or, 
for 7 < uj, the tunneling family. 

It seems likely that most chiral molecules under 
most conditions will be in a regime of strong de- 
coherence 7 > w, see the remarks about D2S2 in 
Scc. lVI Al Whereas it can only be thought of as "chi- 
ral" when in an appropriate rotational state, ammo- 
nia, including its deuterated form ND 3 , is an exam- 
ple of a tunneling molecule in which the transition at 
7/w = 1 can be readily observed in the laboratory. 
Indeed, it appears that it has already been observed; 
see the discussion in Sec. IVIBI One wonders if ad- 
ditional experiments, perhaps using tcchnques other 
than, or in addition to, microwave absorption might 
be helpful in elucidating its behavior near the phase 
transtion. 

In addition to consistent quantum families of his- 
tories we have studied, within the scope of our sim- 
ple model, the flow of information from a tunnel- 
ing molecule to its environment, along with the loss 
of information in the molecule itself. In Sec. [V] we 
used a perspective in which at a later time the infor- 
mation about the quantum state of the molecule at 
an earlier time is thought of as a quantum channel, 
while similar information present at this later time in 
the environment constitutes a complementary chan- 
nel. What happens in both cases depends strongly 
on the type of information considered. Given that 
our model of decoherence, Sec. Ill Bt is based on the 
flow of chiral (X) information — is the molecule left 
or right handed? — to the environment, we were not 
surprised to find this exhibited in our quantitative 
measures, together with a rapid decrease of "com- 
plementary" types of information, corresponding to 
bases mutually unbiased with respect to X, retained 
within the molecule itself. Indeed, there is a direct 
quantitative relationship for short times if one uses 
a Holcvo type of information measure, and an exact 
equality in the instantaneous rates, ([47]) . if in the 
Holevo measure von Neumann entropy is replaced 
with quadratic entropy in order to render the rates 
finite. 

There arc a number of ways in which lines of in- 
vestigation initiated in this paper could be further 
extended. Parity-violation effects could be mod- 
eled as a small energy splitting between chiral states 



[38j | . Also, our model contains only one mechanism 
for decoherence, Sec. [TTJ adding a second would al- 
low a serious comparison with Bcn-Rcuven's formula 
for ammonia as discussed in Sec. IVI Bl Obtaining 
the correct physical interpretation might prove dif- 
ficult given the complexity of the rotational states, 
even for ammonia present as a dilute component in 
a buffer gas. Fluorescence from a two-level atom, 
where decoherence arises from spontaneous decay, 
could be a simpler system for studying the dynami- 
cal phase transition. This transition has been stud- 
ied in terms of correlations among scattered photons 
in |18[ , and it would be of interest to supplement this 



with a description of how the atom itself behaves as 
a function of time. Indeed, even the decay of an iso- 
lated atom initially in an excited state has not, so far 
as we know, been examined using consistent families, 
and studying them might yield valuable physical in- 
sights. 
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Appendix A: Explicit expressions for T and its 
eigenvectors 

The general expression for T(t) = e ts in ([7]) is 



T(t) 



\0 

where, for 7 > uj 



(\ 

e-T*a(t) -e-T*6(i) 
e^ l b{t) -e-T*c(i) 



(Al) 



-2ft 



a(t) = cosh£i + (7/£) sinh£t 

b(t) = (w/£) Binh£i 

c(t) = cosh£rj - (7/O sinh£t 



(A2) 



whereas for 7 < uj 



a(t) = cos jjt + (7/77) sin rjt 

b(t) = (uj/rj) sinrjt 

c(t) = cos rjt — (7/ 77) sin rjt 



r) := \/uP - 



(A3) 



For 7 = uj one has 

a(t) = 1 + yt, b(t) = yt, c(t) = 1 - 7<, (A4) 



16 



where of course 7 could be replaced by w. 

For the eigenvalues Ai and A4 in (fT4]l the left and 
right eigenvectors of T are, trivially, (1,0,0,0) and 
(0,0,0,1), respectively. For A2 and A3, the unnor- 
malized left v and (transposed) right w eigenvectors 
are: 



V2 = (0,£-7, w,0), w 2 = (0,7 
«3 = (0,-7-^,^,0), = (0,7 

where for 7 < u> replace £ with ir\. 
For short times i <?C 1 one has 



T(t) 



for the direct channel and 



e,w,0), (A5) 



/I 











1 









1 - 2 7 i 


\o 










0(i 2 ). (A6) 



T c (t) 



/ 1 o\ 

2/7v^-7 3/2 t 3/2 w/yt 3 / 2 



\l-27i 0/ 

0(t 2 ). (A7) 



for the complementary channel. During this short 
time interval, one can represent the direct channel 
T using using only 2 Kraus operators, instead of 
4 required by the most general qubit channel. The 
physical intuition behind this is that during this time 
interval the system of interest interacts only with 
a qubit environment, being effectively "decoupled" 
from the other environmental qubit (we remind the 
reader that the most general qubit evolution requires 
an interaction with an environment that is repre- 
sented by at least 2 qubits, see e.g. [IH). 



Appendix B: Differential equations for 
consistent families 



As discussed in Sec. IIV A| the forwards and back- 
wards consistency conditions specify how a diameter 
of the Bloch sphere rotates. To determine this for 
the forwards condition, consider the density oper- 
ator which at the initial time is at one end of the 
diameter, and write it in the form 



using the notation of (|4"5|) and (|44|) . Because T is 
unital the master equation (j!2|) for p is equivalent 
to 



dr/dt = S • r, 



(B2) 



or 





" dn~ 










= (ji-S-njn+ 


St 




- dt. 







(B3) 

where S is the lower right 3x3 block of S in ([7|). 
Set r = 1, thus r = n, and take the dot product of 
both sides of (|B3[) with n, noting that n and dn/dt 
are necessarily orthogonal to each other, to obtain: 



dr/dt = n ■ S • n, (B4) 
dn/dt = S • n - (n ■ S ■ n)n. (B5) 



Note that n ■ S • n depends only on the symmetri- 
cal part of S, which is to say the dissipative term, 
proportional to 7, in the master equation Q13p . The 
differential equations in ([30)) are equivalent to (|B5[) 
when n is written in polar coordinates. To obtain 
the differential equations for the backwards consis- 
tency condition, replace S in (|B5[) with S\ corre- 
sponding to the adjoint superoperator T T , and d/dt 
with —d/dt. The resulting differential equations are 
equivalent to (|30|) with 7 replaced with —7, thus 

The consistency conditions are related to the 
motion in the Bloch sphere of the diameter that 
corresponds to the (instaneous) orthonormal basis. 
However, the instantaneous hopping rate k for the 
Markov process can be calculated using the Born 
rule for a very short time interval during which one 
can assume that the diameter remains fixed, as its 
motion (the change in basis) only contributes to 
higher order. When r = 1, k as defined in (|32| 
is equal to (-l/2)dr/dt. Thus, using ([B"4|) . 



K = (-l/2)(n-S-n)=7(l-n 2 ), 



(B6) 



P 



Ul + r-a), r 



(Bl) which, transformed to polar coordinates, is (f3T|) 
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Appendix C: Time-integrated solutions for consistent families 

It is helpful to define fj,(t) := tan <p(t) and v(t) := tan#(i). Then (|30|) and (|33[) can be integrated to give 
the following explicit solutions for the non-stationary consistent families: 

li(t) = m(0) + (w T 2 7 m(0) + ujh(0) 2 ' 



£cosh£i + (±7- w/z(0))sinh£t 



K*) = K0)e^yi± ( 7 /e)^ T ^sinh2Ct + (2 7 /^)( 7T _^||_) silm 2 £i (CI) 

where the top (or bottom) symbol in ± or q= is the solution to (|30|) (or (f33"|0 for the family satisfying the 
forward (or backward) condition. In the case that ui > 7, one can replace every occurrence of £ = -\/7 2 — u> 2 
in (|C1|) with 77 = ^/cj 2 — 7 2 , provided that sinh and cosh are replaced by sin and cos. 



Appendix D: Equality of mutual information and \ measure when forwards consistency conditions 

satisfied 

The key observation is that, when the forward condition is satisfied, T m ,i{P() = J^kQkjPm for each j 
and hence the ensemble of density operators at the channel output commute with each other, so quantum 
(von Neumann) entropies of these density operators become classical (Shannon) entropies in the basis that 
diagonalizes these density operators. Denoting r\ := Ti-(P() and := Tr(P^), we have 

j j 

=S(£ Pj q kj I*/ri) - ^PjSfcQkiPM) 

j,k j k 

j j,k 

=H(££p jqkj r k m /r{} k ) 52 Pj H({ qkj rt/r{} k ) 

j j 
=H(P m ) - H{P m \Pr) = H(P 1 : P m ), (Dl) 

where the k subscript in {•}/. indicates that the set is generated by allowing k to vary. In this derivation, we 
used a property of the von Neumann entropy, for ortho gonal density operators, given on page 513 of |13fl. 
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